BOUNDED HARMONIC FUNCTIONS FOR THE 
HECKMAN OPDAM LAPLACIAN 
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Abstract. We describe the set of bounded harmonic functions for the Heckman- 
Opdam Laplacian, when the multiplicity function is larger than 1/2. We prove 
that this set is a vector space of dimension the cardinality of the Weyl group. 
We give some consequences in terms of the associated hypergeometric func- 
tions. 



1. Introduction 

In this paper we will consider the operator L (called here Heckman-Opdam 
Laplacian) on M" defined, for / a function, by 



(1) = Af{x)+ J2 fc„coth^a„/(a:) 



E 



Here A is the usual Euclidean laplacian, 7?. is a root system, 7?.^ its positive part, 
the Tct's are the orthogonal reflexions associated to the roots and fc is a positive 
function invariant under the action of the r^'s (see the next section). We denote 
by W the Weyl group, i.e. the finite group generated by the r^'s. We denote by L 
the restriction of C to the set of W^-invariant functions. A simpler formula for L is 
given by 

(2) Lfix) = Af{x)+ fc„coth^a„/(x). 

Our main results are the two following: 

Theorem 1.1. Assume that k > 1/2. Then the set of bounded W -invariant har- 
monic Junctions for the Heckman-Opdam Laplacian is exactly the set of constant 
functions. In other words the Poisson boundary of L is trivial. 

Theorem 1.2. Assume that k > 1/2. Then the set of bounded harmonic functions 
for the Heckman-Opdam Laplacian is a vector space of dimension \W\. In other 
words the Poisson boundary of C is W . 

In the next section we will give a precise definition for the terminology " harmonic 
function" . We shall also discuss some consequences of our results in terms of the 
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Heckman-Opdam hypergeometric functions, which are particular eigenfunctions of 
the operator L. 

The first result (Theorem ll.ip was already known for values of fc corresponding to 
the case of symmetric spaces of the noncompact type G/K . The second result (The- 
orem fTT^ is new even for these particular values of fc, but should be also compared 
to the situation on symmetric spaces. There, according to the fundamental work 
of Furstenberg _FJ (see also [GJTp . the Poisson boundary of the Laplace-Beltrami 
operator (but also of a large class of random walks) is K/M . But it was already 
observed that in the Heckman-Opdam (also called trigonometric Dunkl) theory 
the group W often plays the same role than K or K/M . First geometrically, since 
there is a kind of Cartan decomposition: any x € R*^ can be uniquely decomposed 
as w ■ , with the radial part of x (lying in the positive Weyl chamber) and 
w G W. In representation theory also [0 : briefly if H is the graded Hecke alge- 
bra generated by W and the Dunkl-Cherednik operators (see next section), then 
{TC, W) shares some properties of the Gelfand pair (G, K), like the fact that in any 
irreducible flnite-dimensional 7i-module the subspace of M^-invariant vectors is at 
most 1-dimensional. So in some sense Theorem 11.21 is another manifestation (let 
say at an analytical or probabilistic level) of the strong analogy between W and K. 

We should add that the hypothesis /c > is probably sufRcient to get the results 
of Theorem ll.il and 11.21 Here we restrict us to the case k > 1/2, because then the 
stochastic process associated with L (or £) a.s. never hit the walls (the hyperplanes 
orthogonal to the roots, which correspond to the singularities of L), and we need 
it to be sure that the coupling we use is well defined. 

The paper is organized as follows. In the next section we recall all necessary 
definitions. In section [3] we prove Theorem I l.H by using the probabilistic technique 
of mirror coupling. In section [4] we prove Theorem 11.21 by extending the coupling 
to the non-radial process. Our main tool for this is the skew-product representation 
from Chybiryakov [Chyl| , that we have to adapt to our setting. 

Acknowledgments: I warmly thank Marc Arnaudon for having explained to me the 
technique of mirror coupling, and Alano Ancona for enlightening discussions about 
the regularity of harmonic functions. 



Let a be a Euclidean vector space of dimension n, equipped with an inner product 
(•, •), and denote by [} := a-|-ia its complexification. We consider 7^ C a an integral 
root system (see |Bouj ) . We choose a subset of positive roots 7^+. Let = 2a/\a\'^ 
be the coroot associated to a root a and let 



be the corresponding orthogonal reflection. Remember that W denotes the Weyl 
group associated to 7?., i.e. the group generated by the Tq, 's. Let fc : TZ ^ [l/2,+cx3) 
be a multiplicity function, which by definition is VK-invariant. We set 



2. Preliminaries 



ra{x) 



x — {a' ,x)a, 




Let 



a+ {x I Va e 7^+, (a, x) > 0} 
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be the positive Weyl chamber. Let also 0+ be its closure, da+ its boundary and 
Qiog the subset of regular elements in a, i.e. those elements which belong to no 
hyperplane {a = 0}. As recalled in the introduction any a; G a can be uniquely 
decomposed as x = wx^ , with x^ G a+ and w G W. We call x'^ the radial part 
of X and w its angular part. 

For ^ G a, let be the Dunkl-Cherednik operator [O. It is defined, for / G 
C^(a), and x G Oreg, by 

T^f{x)=d^f{x)+ ^ k^-^^^{f(a:)^f{r^x)}~{p,Of{x). 

The Dunkl-Cherednik operators form a commutative family of differential-difference 
operators (see jC] or [O]). The Heckman-Opdam Laplacian C is also given by the 
formula 

n 

1=1 

where {^i, . . . , is any orthonormal basis of a. 

Let A G f). We denote by F\ the unique (see [HO] . [O] ) analytic W-invariant 
function on a, which satisfies the differential equations 

p{T^)Fx — p(X)Fx for all W-invariant polynomials p 

and which is normalized by Fx{0) = 1 (in particular CFx = ((A, A) — \p\'^)Fx)- We 
denote by Gx the unique analytic function on a, which satisfies the differential- 
difference equations (see [O]) 

(3) T^Gx = {X,OGxior all a, 

and which is normalized by Ga(0) — 1. These functions are related by the formula: 

(4) Fxix) = TfJ^ E 

for all X G a and all A G f). 

It was shown in [Sch2j that and are generators of Feller semi- groups 
that we shall denote respectively by (Pt,t > 0) and {P}^ ,t > 0). We will use the 
following definition for harmonic functions: 

Definition 2.1. A bounded or nonnegative function ft, : a — > M is called harmonic 
if it is measurable and satisfies Pth — h for all t > 0. 

Remark 2.1. It is well known that if /i is a function such that Ch — 0, then h is 
harmonic in the sense of Definition 12. II Inversely CoroUarv 12.11 below shows, when 
k > 1/2, that any bounded harmonic function is regular, thus satisfies Ch = 0. On 
the other hand, it is a general fact (which applies for any k > 0), that bounded 
Ty-invariant harmonic functions are regular in a+, but we will not use this fact 
here. 

Observe that by definition Fp is a VK-invariant harmonic function. Moreover it 
is known (see |Sch2| Remark 3.1) that it is bounded. So Theorem 11.11 shows that 
in fact Fp is constant equal to 1. Similarly the functions Gwp^s, for w G W, are 
harmonic and also bounded. This last property follows from Formula since the 
Gwp's are real positive (see |Sch2) Lemma 3.1). In fact one has the following 
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Corollary 2.1. If k > 1/2, then any bounded harmonic function is a linear com- 
bination of the Gwp 's, w W . 

Proof. The only thing to prove is that the G^p's are linearly independent. This 
results from the fact that they are all eigenfunctions of the Dunkl-Cherednik op- 
erators but for different eigenvalues. More precisely, assume that for some real 
numbers {cw)wieWi we have 

^ ^ ^wG^p — 0. 

By applying then the operators p{T(/), with p polynomial, we get 



Cwp{wp)Gwp — for all p. 



^From this, and the fact that Gwp{0) — 1 for all w, it is easily seen that we must 
have = for all w. □ 

3. The W^-invariant case: proof of Theorem 11.11 

In this section we shall prove Theorem 11.11 For this we will use the stochastic 
process [X^ , t>0) associated with L, called radial HO-process, and the so-called 
mirror coupling technique. 

First it is known |Schl| that is a strong solution of the SDE: 

where (-Bt, i > 0) is a Brownian motion on a and 

■= k^a [ coth(a, X^) ds. 

Moreover when k > 1/2, X^ a.s. takes values in a+, or in other words it never 
reaches da+ (see |Schl| ). Now ii x,y G a+, we define the couple {{X^ , Y^),t > 0) 
as follows. Set T = inf{s | X^ = Y^^}. Then by definition X^ is as above, and 
{X^,Y^) is the unique solution of the SDE: 

(5) {Xi^, Y^) = (x, y) + [Bu B',) + {V,\V^\ for t < T, 

where dB^ = rtdBt, with rt the orthogonal reflexion with respect to the hyperplane 
orthogonal to the vector Y^ — X^ (in particular Levy criterion shows that B' is 
a Brownian motion), and 



:= f 



coth(a, YJ^) ds. 



aen+ 

For t > T, we set Y^ — X^ . The existence of this coupling is guaranteed by the 
fact that the SDE ([5]) has locally regular coefficients. We define also by 

._ yiy _ -^w 

and set = \Z^\. It is known [Schlj that a.s. Xf^ /t p, and thus that 
{a,X^) ~ {p,a)t, for all a G TZ'^. From this we see that a.s. sup^o — Vj^l < 
-1-00. Then Tanaka formula ( |RY| p. 222) shows that 

z^ = jt + vt, iort<T, 
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with 7 a one-dimensional Brownian motion and a. s. supj>Q|wt| < +00. In particular 
T is a.s. finite. 

The end of the proof is routine now. Assume that ft, is a bounded W-invariant 
harmonic function. Then it is well known, and not difficult to show, that (h{X^),t > 
0) as well as {h(Y^), t > 0) are bounded martingales. Thus they are a.s. converging 
toward some limiting (random) values, respectively I and I' . Since a.s. = 
for t large enough, we have a.s. / = Then usual properties of bounded martin- 
gales show that 

h{x) =E[;] = ¥\l'] = h{y). 
Since this holds for any x,y S a+, this proves well that h is constant. □ 

4. The non VF-invariant case: proof of Theorem 11.21 

In order to prove Theorem 1 1.21 the first idea is to extend the previous coupling to 
the full process (Xt, t > 0) with semi-group {Pt,t > 0). For this our tool will be the 
skew-product representation founded by Chybiryakov |Chyl| (see jGaYj and |Chy2| 
for the one-dimensional case). Actually Chybiryakov dealt with Dunkl processes, so 
we shall first mention the changes needed to adapt his proof to the present setting, 
and then explain how to combine this representation with the coupling from the 
previous section. 

4.1. Skew-product representation and extension of the coupling. The skew- 
product representation gives a constructive way to define X starting from X^ , by 
adding successively jumps in the direction of the roots. Let us sketch the main 
steps of the construction (for more details see |Chyl| ). First one fixes arbitrarily 
an order for the positive roots: ai, . . . , q;|7j+|. Then for each j S [1, |7?.^|], set 

i<j 

where for any root a, 

U|2 

Decide also that — L. Set 

^'fix) ■.= c-]{x)C'f{x), 

and 

£"+V(x) ■.^c'^lSx)C=f{x). 

The goal is to define inductively a sequence of processes {X^{t),t > 0), j = 
0, . . . , \Ti'^\, associated to the operators C^'s. First X° is just the radial HO-process 
considered in the previous section. Next assume that is the generator of a Markov 
process {X^t),t> 0). Then set 

Jo 

and 

= inf{s > I A^, > t}. 
Using the martingale problem characterization one can see that the radial part of 
X^ is a radial HO-process. Thus for all a G 7?,+ , \{a,Xl)\ > ct, for t large enough 
and c > some constant. In particular the increasing process is bounded. Set 
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= linit^+oo A-j.. Then observe that r/ = +00, when t > . This is essentiahy 
the only difference with the Dunkl case considered in |Chyl| (where was not 
bounded and finite for all t). But one can still see that if 

then killed at time T-', is solution of the martingale problem associated 

with ^J'-J+i (see for instance |EKj exercise 15 p. 263 and section 6 p. 306). The 
next step is to add jumps to X^'^~^^ in the direction of the root aj+i. Namely 
one define a new process X\ also denoted by X^'^^^ ^ ™ |Chyl| section 

2.5, which is solution of the martingale problem associated with Roughly 
X^ is constructed by gluing several paths, all with law or Va^^-^X^'^^^ , such 

that for any two consecutive path the starting point of the second is the image of 
the end point of the first path by the reflexion r-Q^-^j . The lengths of the paths 
are determined by independent exponentially distributed random variables. Here 
the only minor change is that X^ explodes at some time, let say . A change of 
variables shows that 

lim / c-^,(X^(s))ds = +oo. 

t~*T3 Jo 

So for any t >Q, one can define A^ (t) as solution of the equation 

/■*(*) 

c^l^{X^is))ds. 

Jo 

Differentiating this equation one get 

j^A^it)^c^^^AX'{A^t))). 
Then set X^+'^{t) = X^A^{t)), for aU t > 0. The preceding equation gives 

im)^ f c^^^,{x^+\s)) ds, 

Jo 

which in turn shows that X^'^^ is solution of the martingale problem associated 
with C^^^ , as wanted. 

The point now is to combine this construction of X = X^^ ' with the coupling 
of the radial process from section [3l We first take {X''\ with law given by this 
coupling. Then we define the sequence {{X^t),Y^ {t)),t > 0), j = 1,...,|7^+|, 
simply by following the previous construction for the two coordinates. Actually 
this coupling is interesting only when X ~ X^''^ ' and Y — Y^"^ ' never jump, 
but this is precisely what we need. Indeed in this case we have Xt — X'^{t) and 
Yt — Y'^{t), for all i > 0, so they coincide a.s. after some finite time. 

4.2. End of the proof. For any x G a, we denote by the law of {Xt,t > 0) 
starting from x. For e G (0, 1), set 

A^ :— {z a \ Pz[X never jumps] > 1 — e}. 

We know that the process {Xt,t > 0) can jump, so a priori A^ C a. But we know 
also [Schl that a.s. X eventually stops to jump after some finite random time. 
This implies that 

(6) lim P2;[X never jumps after time t] = 1, 

t — *+oo 
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for all x G a. But by using the Markov property, we have for alH > 0, 
Px[X never jumps after time t] — Ej, [PxJ-'^ never jumps]] 
(7) = / P^[X never jumps] dnf{z), 



where /^f is the law of Xt under P^.. So dB]) and ([7]) imply that for all x S a, 
IM^{A^) — > 1, when t — > +cxd. In particular is nonempty. Moreover, by invariance 
of C under we know that for any w G W, the law of {wXt,t > 0) under P^. is P^^;. 
In particular, for any w £ W and any e G (0, 1), we have w{A^ n a+) = A^ Ci wa^. 
Thus all these subsets of A^ are nonempty as well. 

Let now h be some harmonic function. Fix w G W, and take x,y G A^ D wa+. 
Consider the coupling {{Xt,Yt),t > 0) as defined above. Since {h{Xt),t > 0) and 
{h{Yt),t > 0) are bounded martingales, they converge a.s. toward some limits, 
respectively I and I' . We already saw that X^ and a.s. coincide after some 
time. So if both processes X and Y never jump, they must also coincide after some 
time, and in this case we have I = V . Since x,y € A^, this shows that 

\h{x) - h{y)\ = \E[l] - E[l'] \ < 2Ce, 

where C = sup ft.. In particular, by completeness of R, for any sequence {xe)i£(o^i), 
such that Xe & Af n wa+ for all e G (0,1), the limit of h{xe) when e tends to 
exists, and is independent of the chosen sequence. Call 1^ this limit. 

For alH > 0, we denote by Wt the angular part of Xt- Since X eventually stops 
to jump, {wt,t > 0) a.s. converges, i.e. becomes stationary. Then for any w G W, 
define the function hu, on a by 



h^{x)=F, 



lim Wt = w 



By standard properties of Markov processes, we know that these functions are 
measurable, and actually it is not difhcult to see that they are harmonic. Moreover 
the above convergence result for harmonic functions shows that these functions h^, 
w € W, are linearly independent. Then set 



h{x) ^ lyjhy,{x), 



for all a; G a. All that remains to do now is to prove that h — h. Indeed if this 
was true, this would prove that the vector space of bounded harmonic function has 
dimension as wanted. By using the martingale property, we have for any t > 

(8) \h{x) - h{x)\ = \Ex[h{Xt) - h{Xt)]\ < [ \h(z) - h{z)\ dfiUz)- 



We have seen that for all e G (0, 1), 

(9) t^UAe) -> 1 

when t — !■ +00. But it is not difhcult to see (by using the definition of the /u,'s), 
that for any e' > 0, there exists e > such that 

\h{z) - h{z)\ < e' VzeA,. 

Since this holds for any e' > 0, ([8]) and ([9|) show that h = h as wanted. □ 
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Remark 4.1. We have seen in the previous proof that the family {hw)wew is a 
basis of the space of bounded harmonic functions. Since the family {Gwp)wGW is 
another basis, it would be interesting to know the coefficients relating these two 
basis. 
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